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Abstract 

As a candidate for the dark energy, the hessence model has been recently 
introduced. We discuss the critical points of this model in almost gen- 
eral case, that is for arbitrary hessence potential and almost arbitrary 
hessence-background matter interaction. It is shown that in all models, 
there always exist some stable late-time attractors. It is shown that our 
general results coincide with those solutions obtained earlier for special 
cases, but some of them are new. These new solutions have two unique 
characteristics. First the hessence field has finite value in these solu- 
tions and second, their stabilities depend on the second derivative of the 
hessence potential. 

1 Introduction 

In recent years, astronomical observations from type la supernova [1], WMAP 
data [2], and large scale structure surveys [3], have shown that the expansion 
of the universe is accelerated. Although there is no clear understanding of 
the mechanism leading to this acceleration, but it is believed that about 70% 
of the total energy density of universe consists of this unknown energy, i.e. 
dark energy, which leads to this expansion. The simplest explanation of dark 
energy is a cosmological constant A of order (10 -3 ev) 4 . Unfortunately it is 
about 120 orders smaller than the naive expectations, gives rise to the idea of 
a dynamical nature of this energy. The possible dynamical explanations have 
been introduced in different frameworks, such as quintessence [4], phantom [5], 
k-essence [6], tachyons [7], etc. 

In studying the dark energy, the equation of state parameter w,j c = Pde/Pde 
plays an important role, where pde and pde are the pressure and energy density 
of the dark energy, respectively. This parameter is always equal to constant -1 
in cosmological constant model, but it can be a dynamical variable in the above 
mentioned dynamical models. This is an important point since the present data 
seems to slightly favor an evolving dark energy with Wde being below -1 around 
present epoch, [8] , from u>d c > — 1 in the near past [9] . 
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To be definite, we consider the following action 

S = J d^x^—g (--^ + Cde + C rn ) . 1 1 ) 

where g is the determinant of the metric <?p„, 1Z is the Ricci scalar, Cde and 
C m are the Lagrangian densities of the dark energy and matter, respectively. In 
the case of quintessence, the Lagrangian density is 

1 2 

^quintessence = - (dptp) ~ V((j>), (2) 

where is a real scalar field. In a spatially flat Friedmann-Robertson- Walker 
(FRW) universe with homogeneous (f>, w is 

_ 2 /2 - V{<f) 

Quintessence- ^ /2 + y( ^. [6) 

which results — 1 < w qu i n t e ssence- The Lagrangian density of phantom scalar 
field is 

^phantom = ~\ (d^f - V{<f>), (4) 

from which 

-^/2~V(4>) 

phantom " -0V2 + y(0)' [ ) 

where for p > 0, which comes from H 2 = (87rG/3)/?, it results u> p hantom < — 1- 
So we can not cross the phantom-divide- line w — — 1 in quintessence or phantom 
model alone. A possible way to overcome this problem is considering two real 
fields, which one behaves as quintessence and other one as phantom field. The 
resulting model, called the quintom model, has the following Lagrangian [9, 10] 

1 2 1 2 

Ajuintom = g (<9 M 0l) -^(0^2) -V(<t>l,i>2), (6) 

with 

_ ^i 2 /2-fo 2 /2-y(fc,fe) 

^quintom — .9 .9 ' v / 

0! /2-(t> 2 /2 + V(<H,<h) 

• 2 -2 

Now it is obvious that w qu intom > —1 when (pi > (f>2 and u> qu intom < — 1 

• 2-2 

when <fii < (f>2 ■ So crossing the phantom-divide-line is, in principle, possible 
in quintom model. See, for example, [10] and [11]. 

Instead of introducing two independent real scalar field to describe a quintom 
model, it is also natural to consider a single complex scalar field. The resulting 
spintessence model of dark energy [12-16], has the following Lagrangian density 

£spintessencc = \ (0"$*) - U(|$|). (8) 

Using $ = 0i + i(p2, the kinetic term of eq.JHJ reduces to the kinetic terms of 
eq.lJSJ. Also the above Lagrangian is invariant under $ — > e la & which leads to a 
conserved charge. Unfortunately this model suffers from the problem of Q-ball 
formation [12,14] . Q-ball is a kind of nontopological soliton which except in some 
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special cases of spintessence with unnatural potentials, grows exponentially and 
depends on the potential, can be either stable at the late-time to be a dark 
matter, or decay into other particles. Therefore the spintessence model can not 
be a viable candidate for the dark energy. 

To avoid the difficulty of Q-ball formation and also to introduce another 
possibility for mysterious dark energy problem, a non-canonical complex scalar 
field, called hessence, has been recently introduced in [17]. In the hessence 
model, the phantom-like role is played by the so called internal motion 9, where 
is the internal degree of freedom of hessence. There is a conserved charge Q 
in this model which makes the physics of hessence more interesting, and the 
transition from Wh > — 1 to Wh < — 1 or vice versa is also possible. Another 
interesting feature of hessence model is that it is free of big-rip [18]. If w < — 1 
in an expanding FRW universe, then the positive energy density of a phantom 
matter generally becomes infinite in finite time, overcoming all other forms of 
matter and hence leads to the late-time singularity called the "big-rip" [19]. 

By considering two specific hessence potentials, i.e. the exponential and the 
(inverse) power law, and four different interaction forms between hessence and 
background perfect fluid, the late-time attractors of hessence model have been 
studied in [18]. In each case, different scaling and hessence-dominated solutions 
have been obtained and their stability properties have been studied. 

In this paper we are going to study the late-time attractors of the almost 
general hessence model, with arbitrary hessence potential and almost arbitrary 
hessence-background matter interaction term. By almost arbitrary, we mean 
that the hessence potential and hessence-background matter interaction terms 
are arbitrary functions of dimensionless variables defined in (|22(l . Specifically, 
we mean eas. (!33l) and (1341) . We show that there always exist some stable so- 
lutions, in scaling or hessence-dominated form, which some of them have not 
been appeared in special cases studied in [18]. These new solutions have two 
interesting unique characteristics which are absent in other solutions. First their 
stability depends on the second derivative of the hessence potential, and second 
the hessence field has finite value in these solutions. The significance of the sec- 
ond derivative of the potential in the late-time behaviors has been also revealed 
for quintessence model in [20] , in which some conditions have been imposed on 
the first and second derivatives of the potential. 

The scheme of the paper is as follows. In section 2, we briefly introduce the 
main points of hessence model and the system of equations which determines 
the critical points in terms of dimensionless variables. In section 3, we consider 
the general hessence potential, but assuming no hessence-background matter 
interaction. It is shown that there exists five general solutions for critical points 
which three of them are stable under specific conditions. Finally in section 4, the 
hessence potential and hessence-background matter interaction are considered 
arbitrary and it is shown that there are, in general, six classes of solutions for 
critical points. The stability of solutions is discussed in special cases. It is shown 
that all the solutions of [18] can be obtained from our general results. 

We use the units h = c = 1, k 2 = 8irG and adopt the metric convention as 
(+, — , — , — ) throughout the paper. 
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2 Hessence model 

Following [17], the hessence field introduced by a non-canonical complex scalar 
field 

* = 0i (9) 

with Lagrangian density 

c h = \[ {d^f + (d^*) 2 } - u(<p 2 + $* 2 ) = \ [ {d^f - tf(d,ef ] - 

(10) 

where the new fields (</>, 9) are defined through 

4>i = </>cosh#, <f>2 = (j)smh6. (11) 

In a spatially flat FRW universe with scale factor a(t), the equations of motion 
for <f) and 9, when they are considered homogeneous, are 

4> + 3Hj> + 4>e 2 + = o, (12) 

4> 2 9 + (24><p + 3H(/) 2 )9 = 0, (13) 

where H = a/a is the Hubble parameter and overdot and subscript " ^" denote 
the derivatives with respect to cosmic time t and <f>, respectively. Eq. (|13|l implies 

Q = a 3 <f) 2 9 = const., (14) 

where Q is the total conserved charge due to the symmetry of Lagrangian 11U|) 
under the transformation <fi —> <fi and 9 — ► 9 — ia. Substituting ca. (|14fl into (|12fl . 
one has 

a 6 (j> 3 

The pressure and energy density of hessence are 



The Friedmann equation and Raychaudhuri equation are given by, respectively, 

k 2 

K 2 

# = — 2^ Ph + Pm + Ph + Pm ^ ' 

where p m and p m are the pressure and energy density of background mat- 
ter, respectively. The background matter is described by a perfect fluid with 
barotropic equation of state 

Pm = W m Pm = (j ~ l)p m , (19) 

where < 7 < 2. In particular, 7 = 1 and 4/3 correspond to dust matter and 
radiation, respectively. 
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To introduce the interaction between hessence and background matter, it 
is assumed that it can be described by an interaction term C in the energy 
balance [18,21] 

p h +3H{p h +p h ) = -C, (20) 
f J rn m ~T~ Pm ) = C, (21) 

which preserves the total energy conservation equation p to t + 3H (ptot + Ptot) = 
0. C — corresponds to no interaction between hessence and background matter 
and when C ^ 0, a new term due to C will appear in the right hand side of 
en.lfTKjl, 

Following [22] and many other papers, if we introduce the following dimen- 
sionless variables 

K(j) K,VV K\/~Pm K Q 

V = ~ 7S>7' z = /77„ i u = — T> 17 = ^77^77 > ( 22 ) 



V6i?' V3i?' ' ' %/3# ' " k0' VGHa 3 

then using eas. (|16|l - l|18|l . (|20p . and J2U, the evolution equations of these variables 
become 

x 1 = 3x (x 2 -v 2 + |z 2 - l) - uv 2 - y|y 2 / - C 1} (23) 

y' = 3y (x 2 - v 2 + |z 2 ) + \J\xyf, (24) 

z' = 3z (x 2 -v 2 + |z 2 - |) + C 2> (25) 
v! = -xu 2 , (26) 
v' = 3v [x 2 -v 2 + ^z 2 - l) - (27) 

Prime denotes derivative with respect to the e-folding time TV = In a, and 

/=% (28) 

d EE -^L, C 2 EE — ^ = -d. (29) 

The Friedmann equation (|17fl becomes 

x 2 + y 2 + z 2 - v 2 = 1, (30) 
and the fractional energy densities are 

r> P h 2,22 r> ^ m 2 / i\ 

S2/i = — =x+y-w, il m = — = z, (31) 

Pc Pc 

where p c = is the critical energy density. The equation of states of hessence 
and whole system are 

Ph X 2 -V 2 -y 2 Ph+Pm 2 2 2,/ i\2 /o A 

p h x 1 - v 2 + y z p h + p m 

The critical points (x, y, z, u, v) are obtained by imposing the conditions x' = 

y' = Z f = U' = V' = 0. 
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Note that we take 



/ = f(u), 



(33) 



and 

C\ = Ci(x, z, u), C 2 — C 2 {x, z,u). (34) 

It is because V = V(4>), so / is assumed to be a function of only one variable u. 
Also because of ea. ((30(l . only four of the variables ((22(1 are independent, which 
we can take them x, y 7 z, and u. But the dependence of hessence-background 
matter interaction to the potential V (or variable y) is meaningless, so Cs are 
taken to be arbitrary functions of variables x, z and u. Really ea . 1)33(1 does not 
constraint the potential V, but ea . ((34(1 restricts the possible interaction term C. 
In this way the eas. ((23(1 - 1(271) become autonomous and we need not consider any 
further variables. For example / is not and extra variable since /' = leads to, 
for an arbitrary potential, u' = 0. 

3 Attractors in C = case 

To obtain the attractors for arbitrary hessence potential and when there is no 
hessence-background matter interaction, we must solve ea. 1(30(1 and the set of 
equations (|2^|l - (|2T)l . when setting zero, in C\ = C 2 = 0. Ea. (|2l)|) results u = 
or x — and ea. l(2"5|) results z — or x 2 — v 2 + (j/2)z 2 — 7/2 = 0. So we have 
four possibilities: 1= (u = 0,z = 0), 11= (u = 0, x 2 - v 2 + (7/2) z 2 -7/2 = 0), 
111= (x = 0,z = 0) and IV= (x = 0,x 2 - v 2 + (~f/2)z 2 - 7/2 = 0). In type I 
solution, ea. (|2Tjl . using l(3^|) . reduces to vy 2 = 0, so it divides to 1.1= {u = 0, z = 
0, v = 0) and 1.2= (u — 0, z — 0, y — 0) solutions. The remaining variables can 
be easily found. The final results are represented in Tabled 



Label 


Critical Point (x, y, z, u, v) 






w h 


W eff 


I.l 

1.2 


-J- v / 1 — £-ooo 

x 2 > 1, 0, 0, 0, ±Vx 2 - 1 


1 
1 






f 2 

-1+*- 
T 3 

1 


T 3 
1 


II 

III 

IV 


0, 1, 0, / = o, 

0, 0, 1, any, 


3j 
P 
1 




1 _ 52 
1 p 


1 


-1 + 7 
-1 
any 


-1+7 

-1 
-1+7 



Table 1: Critical points for arbitrary hessence potential when there is no 
hessence-background matter interaction. 

The solutions with flh — 1, i.e. solutions I.l, 1.2 and, III, are hessence- 
dominated, solution IV is background-matter-dominated and solution II is scal- 
ing solution. In solution III, u must be found by solving / = f(u) = 0. Note 
that in all cases, w e // > — 1. It is also interesting that for potentials where 
f 2 {u = 0) = 37, the solutions I.l and II become degenerate. 

In examples considered in [18], the potentials are 

V x = V e- X ^, (35) 
V 2 = Vb («</>)", (36) 
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or, using (|2*81l . 



h = -\ 

n nu 



(37) 
(38) 



It can be easily checked that our solutions (1.1, 1.2, II and IV) and (1.1, 1.2, and 
IV) reduce to those obtained in [18] for V\ and V2, respectively. The solution 
II does not exist for potential V2 since /^(u = 0) = 0. The solution III is a 
new solution which has not been appeared in [18]. This is because the equation 
/ = results A = for Vi, which is not acceptable, and results u = for V2, 
which reduces solution III to 1.1. 

To study the stability of the critical points I.l-IV, we must consider a small 
perturbation about the critical point (x,y, z, u, v): x — > x + 5x, y — > y + 8y, 
z —* z + Sz, u — > u + Su, and v — > v + Sv, in eqs. J23J-(EHJ w hh C\ = C2 = 0, 
which due to Fricdmann constraint 13UII . only four of them are independent. In 
this way one can found a 4 x 4 matrix M defined through 



d 
<W 



( Sqi \ 



Sq 1 
Sq2 

V ^94 



= M 



J 



Sqi 
5q2 
Sq3 

V ^94 



(39) 



/ 



where (q%, ■ ■ ■ , q^) are four chosen independent variables. The critical solutions 
are stable if the real part of all the eigenvalues of matrix M are negative. The 
eigenvalues of matrix M for our solutions are as following: 



1.1 : (-3y 2 ,-6y 2 ,i(/ 2 -3 7 ),0), 
I.2:(0,0,|(2-7),3 + i/|x/), 



3 ^9 + 2V6m 2 /' 3 ^9 + 2V6u 2 /' 
_ 2 + 2 ,_ 2 2^ 

IV: (o,^,3(7-2)^( 7 -2) 



(40) 
(41) 

), (42) 
(43) 



in which ■ ■ ■ , 54) = (a;, y, z, u) and /' = (d//du) fi . It is clear that 1.1 solution 
is stable if 

P < 37, (44) 



1.2 solution is unstable since 7 < 2, III solution is stable if 



du 



< 0, 



(45) 



and IV is an unstable solution since 7 > 0. 

For solution II, it is easier to use (q±, • • • , 34) = (x, z, u, v). Then it can be 



easily found that l\ = 0, I2 
I 2 + bl + c = Q with 



1(7 — 2), and ^3 and Z4 are roots of equation 



5(7-2), 
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c = r_(2- 7 )(.f 2 -3 7 ). 



(46) 



Is and I4 are non-positive if b > and c > 0. As 7 < 2, the solution II is stable 



So for any potential V , there always exists at least one stable attractor. A 
hessence-dominated attractor (1.1) if f 2 < 37 or a scaling attractor (II) if f 2 > 



It is interesting to note that the solution III has two unique properties. First, 
it is the only stable attractor which has the non-vanishing u value, i.e. finite 
value of hessence field (f>. Second, it is the only attractor which its stability 
depends on the derivative of / (eq.(35J). In other words, the stable attractors 
1.1 and II can not distinguish between different potentials with the same / value, 
but the attractor III does. 

As examples of potentials which have the solution III as a stable attractor, 
we may consider V3 = Vo sin(K</>/v / 6) and V4 = Vo cos(k0/v / 6). For V3, we have 
fa = — (1/ku 2 ) cot(l/w), which results u — [(2n+ l)7r/2] _1 , with n = 0, 1, 2, • • • , 
as the solution of fa — equation. It really has infinite number of attractors. 
Then (d//du) a = — (l//c)[(2n+ l)7r/2] 4 < 0, which shows that the attractors 
are stable. The same is true for V4 potential. 

4 Attractors in the presence of hessence-background 
matter interaction 

In this case, we must solve ea. (|3U|l and the set of equations (|23 p - H27[l . when 
setting zero, for arbitrary C function. Eq. l|26|l results u = or x — 0. In each 
case, we consider eight cases in which each of the variables y, z and v has two 
possibilities, zero and not zero, and then check the consistency of the equations. 
The final results are as following: 



solution l:{x= — - r , y=\ l- J —, z = u = v = 0, d = C 2 = 0}. (48) 



Note that the equations C\ = C2 = imply that the functional form of C\ and 
C2 must be such that they are identically equal to zero at this critical point, 
otherwise this solution does not exist. 

solution 2 : {x 2 > 1, y = z = u = 0, v = ±\/x 2 - 1, C x = C 2 = 0}. (49) 

In this case, the equations C\ = C2 = can generally determine the allowed 
value of x. If C\ and C2 are identically equal to zero at y = z — u = 0, as they 
are in C — case, then x can choose any arbitrary value. 



if 



P > 3 7 . 



(47) 



3 7 . 




solution 3 : {x — y = 0, z 



1, u, v = 0, d = C 2 = 0}. 



(50) 



The value of u is generally determined by solving C\ — C'2 — 0. 



solution 4 : {x = 0, y = 1, z — 0, u, v — 0, C\ 




If, C 2 - 0}, 



(51) 
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where the last two equations can generally determine u. 

„_,_ 2 ,7-2 7 

= —3z(x+ — z 

V 2 2 

(52) 



solution 5 : {y = u = v = 0, d = 3x{x 2 + ^z 2 -l), C 2 = ~3z(x 2 + ^z 2 -^), x 2 +z 2 = 1}. 



x and z are found by solving the above equations. The last solution is: 



solution 6 : { y = \Jl - x 2 - z 2 , u = v = 0, C x = 3x(x 2 + - 1) - J^ff, 



=2 , 7=2 7x Q ,-2 , 7=2> , / 3 ; 



C 2 = -32(^ + ^ - 3(* 2 + ^) + y -5/ = 0}. (53) 

At C = 0, solutions 1, 2, 3, 4, and 6 reduce to 1.1, 1.2, IV, III, and II of 
Tabled respectively. The solutions 1, 2, and 4 are hessence-dominated, solution 
3 is background-matter-dominated and 5 and 6 are generally scaling solutions. 
Among these solutions, there are only two solutions 3 and 4 in which u can 
principally be different from zero, which have not been appeared in examples 
discussed in [18]. 

In [18], besides the non-interacting C — case, three following interactions 
have been considered: 

-o(ii) /3 2 

1 = V 2 aZ ' 
n (in) _ 3/3 

1 "2? 

Cf v) = (54) 

in which a, /3 and rj are some constants. In both of the solutions 3 and 4, x 
is zero, so C^ 111 -* diverges and C^* is not generally well-defined and therefore 
these critical points do not exist in these cases. So we only consider the C± 
case. For V — V\ potential, with f\ = —I, solution 3 does not exist since 
C\ = leads to a — which is not acceptable, and solution 4 does also not 

exist as C\ — —\j\f results in I = which again is not acceptable. For V = V% 

potential, with fi = nu/yE, solution 3 leads to a = which is not acceptable, 
and solution 4 results in u — 0, which does not lead to a u ^ solution. 

It may be useful to reproduce all the critical points of at least one of the 
cases studied in [18] in more detail. We consider, as an example, V = V\ and 
C = C[ U) . So 

/ = -I, Ci = \j^az 2 , C 2 = J^axz. (55) 
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The solutions 1-6 result in: 



I / I 2 

solution 1 : {x = — j=, y = y 1 — — , z — u — v — 0}, 



solution 2 : {x 2 > 1, y = z = u = 0, i> = ±y x 2 — 1}, 
solution 3 : C\ — — > a — (is not acceptable), 

/3 - 

solution 4 : Ci = — y — / — ► i = (is not acceptable), 



/2 a / 2a 2 



, . r _ /3 7 /2a 2 -3(7-2) 7 + 2aZ _ 1(1 + a) - 3j 

solution 6 : {a; = i / — ,y=\ —. ^5 i z =\ — r, ^ — > 

V 21 + a V 2(a + l) 2 V (Z + a) 2 



which coincide with those in Table II of [18]. Note that in that table, four of 
the solutions ( 2p, 2m, 4 and 5) are not independent solutions and are special 
cases of the first solution. 

The stability studies of these critical points depends on the precise value of 
the function C . But it may be interesting to study the conditions under which 
the derivative of / becomes important in the stability properties of the critical 
points. Consider the most general case C\ = C\(x, z,u). It can be shown that 
the coefficient of /' term in equation det(M — 11) = 0, where 1 stands for 4x4 
unit matrix, is 

~ 2 u 2 

-^-h{x,y,z,u), (62) 
z 2 - 

which shows that the derivative of / survives only if 

y ^ and u ^ 0. (63) 

Therefore the solution 4 is the only one in which we expect /' may be important. 
Focusing on this solution, the eigenvalue problem is still very complicated in 
general case and it is better to consider the more restricted cases. As an example, 
we consider the potentials and C's in which / = [C\j z) = (dCi/dx) — 0. Under 
these conditions, the equation det(M — 11) — leads to 

I 4 + C1J 3 + a 2 l 2 +a 3 l + a 4 = 0, (64) 

where 

3 

ai = 9 + -7, 

27 - A3 - 

aa = 18 + _ 7 -u2(C liT , + y-/'), 

a 3 = -3u 2 (<?!,„ + y|/')(2 + \) + 27 7 , 

a4 = -97" 2 (Ci,« + \/!/'), (65) 



where C\. u = (dCi/du). Now it is well known that li = — ai, Ylj-^j hlj = &i, 
^2i=£j^k liljlk — ~ a 3i an( i hhhh — a 4j where l\ ■ ■ ■ I4 are the roots of ea.l]64p. 
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So all the roots are non-positive only if <ij > {i = 1, • • • , 4), which provided if 



In C = 0, this condition of stability reduces to ea. (|4*5*)) . 

5 Conclusion 

In this paper, we study the attractor solutions of the general hessence model by 
studying the 4-dimensional phase space of the theory. The hessence model is a 
non-canonical complex scalar theory which can be a candidate of dark energy 
with some interesting properties, among them is the possibility of crossing the 
w = — 1 line. Comparing the Lagrangian of hessence model with quintom model 
shows that the 0-( or equally charge Q-)term plays the role of the phantom field. 
In Q = 0, the model reduced to quintessence model with no w = —1 crossing. 

We consider an arbitrary hessence potential V(4>) and almost arbitrary 
hessence-background matter interaction term C, and find several results. We 
show that in C = 0, there is always at least one stable attractor which depends 
on the value / = f(u = 0): For f 2 < 3j the hessence-dominated attractor 1.1 
and for f 2 > 37 the scaling attractor II of Table 1. In all the attractor solutions 
we have (f> — > 00, except for the attractor III where cf> is finite and also its stabil- 
ity depends on the derivative (df/du)u- This kind of attractor did not appear 
in the previously studied cases and can be seen in potentials like sin(K</>). v is 
zero in all the stable attractors of Table 1, but for arbitrary C-term, this is not 
the case. 

For general C, we show that there generally exist six classes of attractor, 
which all of them can be stable in special cases ( for example V = Vi and 
Ci = C'^ 1 ^ of [18]). Among these solutions, solution 2 is the only one with 
property v =/= 0. The solutions 3 and 4 have finite <f> value and the the solution 
4 is the only one which its stability depends on (df /du)u, the features that can 
not be seen in the previously studied potentials. 

Acknowledgement: We would like to thank the research council of the 
University of Tehran for partial financial support. 
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